Can quarkonia survive deconfinement? 
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We study quarkonium correlators and spectral functions at zero and finite temperature in QCD 
with only heavy quarks using potential models combined with perturbative QCD. First, we show that 
this approach can describe the quarkonium correlation function at zero temperature. Using a class of 
screened potentials based on lattice calculations of the static quark-antiquark free energy we calculate 
spectral functions at finite temperature. We find that all quarkonium states, with the exception 
of the IS bottomonium, dissolve in the deconfined phase at temperatures smaller than l.STc, in 
contradiction with the conclusions of recent studies. Despite this the temperature dependence of 
the quarkonium correlation functions calculated on the lattice is well reproduced in our model. 
We also find that even in the absence of resonances the spectral function at high temperatures is 
significantly enhanced over the spectral function corresponding to free quark antiquark propagation. 

PACS numbers: ll.15.Ha, 12.38.Aw 



I. INTRODUCTION 

The study of quarkonia at finite temperature is in- 
teresting for several reasons. First, due to their small 
size, these heavy quark-antiquark bound states provide a 
bridge between perturbative and nonperturbative Quan- 
tum Chromodynamics (QCD), testing forces at interme- 
diate distances. At high temperatures, color screening, 
which is usually understood in terms of in-medium mod- 
ification of interquark forces, occurs. Based on this Mat- 
sui and Satz argued that above the transition temper- 
ature Tc screening efi'ects are strong enough to lead to 
dissolution of the J/tp state. This can then be used as 
a signal of quark gluon plasma formation in heavy ion 
collisions [l|. 

Because of the large quark mass m = rricb 3> Aqcd, 
the velocity v of heavy quarks in the bound state is small, 
and the binding effects in quarkonia at zero temperature 
can be understood in terms of a nonrelativistic potential 
model with a Coulomb plus linear form, known as the 
Cornell potential [2, i3h A: lS] • Potential models appeared 
to be very successful in describing the quarkonium spec- 
trum and have been extensively used in the past 20 years, 
see Ref. [6]. More recently, an understanding has devel- 
oped on how to derive the potential models from QCD 
using a sequence of effective field theories : nonrelativis- 
tic QCD (NRQCD), an effective theory where all modes 
above the scale m are integrated out, and potential non- 
relativistic QCD (pNRQCD), an effective theory in which 
all modes above the scale mv are integrated out 0, Q ■ 
The concept of the potential can be given a solid field 
theoretical definition in this framework at any order of 
perturbation theory. In particular, the quark antiquark 
potential is defined in terms of the expectation value of 
the Wilson loops. Rclativistic corrections to the poten- 
tial can also be calculated. They are expressed in terms 
of the Wilson loops with appropriate insertions of electric 



and magnetic fields, see Refs. [ol. [lOl. . 

Based on the success of the potential model at zero 
temperature, and the idea that color screening implies 
modification of the interquark forces, attempts to under- 
stand quarkonium properties at finite temperature using 
potential models have been made [U, [H, [13, [3| ■ While in 
these works phenomenological potentials have been used, 
more recent studies went one step further and attempted 
to connect the potential to lattice calculations of the fi- 
nite temperature free energy of a static quark antiquark 

pair [n, til. Hi [n, m ioiii, m . 

Quarkonium states at zero temperature are well de- 
fined: their widths are very small compared to their bind- 
ing energies, at least for states below the threshold. At 
finite temperature the situation is different. We expect 
that all quarkonium states will acquire a sizable width, 
which increases with increasing temperature. At some 
temperature the width becomes large enough that it is no 
longer meaningful to talk about individual quarkonium 
states. Instead, one should consider the spectral func- 
tion, which contains contributions from all possible states 
in a channel with given quantum numbers. Furthermore, 
the spectral function in the vector channel is a quantity 
which can be measured directly, since it is proportional 
to the dilepton production rate. Quarkonium spectral 
functions at finite temperature have been considered only 
relatively recently. Using lattice QCD, charmonium cor- 
relators have been calculated and the corresponding spec- 
tral functions have been extracted u sing the Maximum 
Entropy Method (MEM) [H, Hi, HE Hi 113, HI US 113 • 
Although this approach can in principle provide the ulti- 
mate solution to the problem of in-medium quarkonium 
properties, current calculations have serious limitations 
and cannot give detailed information about quarkonium 
spectral functions. Using the MEM at zero tempera- 
ture one can reconstruct the basic features of the spectral 
fimctions: the ground state, the excited states, and the 
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continuum [30|. Individual excited states, however, can- 
not be resolved. At finite temperature, even resolving 
the ground state appears to be difficult with existing lat- 
tice data. The only statement that can be made at this 
time, is that in the pseudoscalar channel, quarkonium 
spectral functions do not show significant changes up to 
temperatures as high as l.STc, while in the scalar channel 
the spectral function is strongly modified just above the 
transition temperature [30| . 

For this reason, quarkonium correlation functions have 
been studied using a simplified model of the spectral 
function, which contained bound states and a perturba- 
tive continuum [sil, [13] . The results of these calculations 
have been compared to lattice results and no agreement 
has been found. Very recently this approach has been 
extended using the full nonrelativistic Green's function 
of heavy quark antiquark pairs to estimate the spec- 
tral function and the corresponding Euclidean correla- 
tors [21I, I23, HI]. However, even in these approaches no 
agreement with lattice calculations has been found. 

In the present paper we study spectral functions in 
the pseudoscalar, vector, scalar and axial-vector chan- 
nels which correspond to ric {rjb), J/ip (T), XcO (Xbo) 
and Xci ixbi) charmonium (bottomonium) states, respec- 
tively. In the energy region below and near the contin- 
uum threshold the spectral function is calculated using 
a potential model and nonrelativistic Green's function. 
Since the potential at finite temperature is not known, 
we consider a class of screened potentials based on lattice 
results on the static quark-antiquark free energy. Well 
above the threshold, the nonrelativistic spectral function 
is matched to the perturbative fully relativistic result. 
From this the Euclidean correlators are calculated. We 
then compare these correlators to the results of recent lat- 
tice calculations. We find that the lattice data does not 
necessarily imply survival of different quarkonium states. 
Rather, despite the fact that most quarkonia, with the 
exception of the IS bottomonium, are dissolved at tem- 
peratures smaller than l.STc, agreement with the corre- 
lator data is found. Clearly, dissociation temperatures 
previously quoted in the literature (e.g. 2Tc for J/ijj ) 
have not been seen in our analysis. 

The rest of the paper is organized as follows: The 
framework for calculating quarkonium spectral functions 
is discussed in Section II. In Section III we show our anal- 
ysis of the quarkonium spectral functions and Euclidean 
correlators at zero temperature and compare to lattice 
QCD results. Sections IV and V contain our results on 
the finite temperature spectral function and correlators. 
Finally, in Section VI we give our conclusions and out- 
look. The reader not interested in technical details can 
skip Sections II and III. Further technical details of our 
calculations are presented in the Appendices. 



II. QUARKONIUM SPECTRAL FUNCTIONS IN 
THE NONRELATIVISTIC LIMIT 

In this section we discuss quarkonium spectral func- 
tions in the free theory as well as in the nonrelativistic 
limit. In what follows we consider the case of zero spa- 
tial momentum, i.e. quarkonium at rest. The spectral 
function is defined as the imaginary part of the retarded 
current-current correlator 

a{uj) =--lmDR{uj), (1) 

TT 

zDniio)^ J dte'^'e{t) J d'x{[j{t,x),j{0,0)]), (2) 

It carries information about all the possible states with 
a given quantum number, which is fixed by the current 

j = qilq (3) 

with fl — 1,75,7^,7^75 for scalar, pseudoscalar, vector 
and axial vector channels respectively. 

For large uj the spectral function can be calculated in 
perturbation theory (see e.g. [sl-lssj) 

Here a = 1, 6 = for pseudoscalar, a = 2, 6 = 3 for 
vector, a — l^b = —1 for the scalar and a = 2,b — 
—2 in the axial- vector channel, respectively. In leading 
order perturbation theory the threshold is sq — 2mc.b- 
The coefficient C of the leading perturbative correction 
has been calculated only for the massless case [36]. The 
number of colors in QCD is iVc = 3. 

While perturbation theory is reliable away from the 
threshold, a; ^ So, the physics becomes quite com- 
plicated near the threshold, even in the weak coupling 
regime [13, [H, [H, [^] ■ Close to the threshold the quark 
and antiquark move slowly allowing enough time for mul- 
tiple gluon exchange. Adding an extra gluon exchange 
does not lead to a suppression by as- In this case we 
need to resum ladder diagrams. In the following we dis- 
cuss this resummation separately for the pseudoscalar 
and scalar channels. The vector channel has been dis- 
cussed in detail in Ref. [!§] , while the axial- vector chan- 
nel is completely analogous to the scalar case. 



A. pseudoscalar channel 

The summation of ladder diagrams corresponds to 
solving the integral equation for the vertex function 

(5) 
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and inserting the solution of this equation into the 1-loop 
expression of the meson correlator (see e.g. [39]) 



Nr 



(2^ 



Tr 



Sf{p + |)r(p, q)SF{p ~ |)75 



(6) 

In general this task is complicated, but it simplifies con- 
siderably in the nonrelativistic limit. In this case follow- 
ing Ref. [3^ we can replace the quark propagators with 
the corresponding nonrelativistic forms 



Sf (p 



Si 



(- 



(1 + 7o)to + 70 (po + E/2) + 7 ■ p 
2m{§+po-i, + ^e/2) ' 

(1 - 7o)m + 7o(po - E/2) + 'j ■ p 
Ml -Po-^ + ^e/2) 



(7) 



Here we have taken into account that q — {2m+E, 0) and 
E <^ 2m. The 1-gluon exchange operator in this limit is 

V{k^p) = V{\k-p\) - |a,4^I?oo(|^-pl), where i?oo is 
the temporal Coulomb gauge gluon propagator. In this 
limit the vertex function can be chosen to be independent 
ofpo M- 

The leading order result for the pseudoscalar spectral 
function can be obtained by replacing F = 75 in Eq. Q 
and using the nonrelativistic form of the quark propaga- 
tors above with only the first term in the denominator. 
The retarded nature of the meson correlator D{q^) is en- 
sured by the +ie prescription in the quark propagators. 
This gives for the free nonrelativistic spectral function 

a{E)=--lmD{q^) = ^m^/^E^/\ (8) 

Note, that this result can also be obtained from Eq. (|4]) 
by writing lo = 2m + E and sq = 2m, and expanding in 
E /m to leading order. 

Defining the scalar function 



(1 + 70)^.. (1 - 7o) 
r(p, E) 



7y 75 7y r(P, E) 

(9) 

and performing the integral over fco explicitly, the integral 
equation for the vertex function can be written as 



f{p,E) = l 



1 



(27r)3£; + ie-fc2/r 



-V(\p-k\)t{k,E). 



By introducing the nonrelativistic Green's function 

1 



(10) 



G"''ik,E 



E 



T{k,E), (11) 



the equation for the scalar vertex function r(fc, E) can 
be rewritten in the form 



E 



P_ 

m 



G"'^(p,S 



d^k 
(2^ 



Vi\p-k\)G'^^ik,E + ie), (12) 



which is the Schrodinger equation in momentum space. 
In a more familiar form, in coordinate space, it reads 



-— + T/(r) - (E + ie) 
m 



Here 

G"''(f,P,£: + ie) 



G'^''{r,r',E+ie)^5^{r-r'). 

(13) 



d^k 



i2n) 



3e~*'=-(*'^''')G"''(fc,£: + ie) 



(14) 

is the nonrelativistic Green's function in coordinate space 
and 



Vir) 



d^k 



ik-r 



(2^)3 

is the potential. Thus we can write 



(15) 



a{E) = ^Im 

TT 



d^k 

j^G-ik,E^^.) 



2Nr 



ImG"'-(f,r',i?-|-ze)|^_p^„. 



(16) 



Therefore, in order to calculate the pseudoscalar spectral 
function in the nonrelativistic limit we have to solve the 
Schrodinger equation (fT3| for G"'''(r, r',E + ie) and take 
the limit f = r' = 0, according to ((T| 



B. Scalar Channel 

The calculation of the scalar spectral function is some- 
what more complicated. To understand the problem bet- 
ter, let us consider the noninteracting case first. From the 
structure of the quark propagators in the nonrelativistic 
limit it is clear that for the scalar vertex 17 = 1 the meson 
correlator is equal to zero in leading order of a 1/m ex- 
pansion (c.f. Eqs. ([7])). Therefore the second and third 
terms in the numerator of Eqs. ([7]) should be retained. 
We then get 



D{q' 



Nr. 



d^P 



^-E/r 



{2-)^(f +Po + .f-^)(|+Po+^|-^) 



(17) 



Taking the imaginary part and performing the integral 
we get 



a{E) 



27r2' 



^1/2^3/2 



(18) 



Note that we arrive at the same result if we consider 



the nonrelativistic scalar vertex f2 



1-P 



instead if the 



relativistic one = 1. In fact, this type of nonrelativistic 
vertex is used to study the Xc,b states in NRQCD [4l| [70j . 
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Repeating all the steps discussed in the previous sec- 
tion we write the correlator 



Tr 



SFip+^)T{p, q)SF{p- |)^ 

(19) 



where the vertex function T(p, q) satisfies the equation 



r(p,9) 



1-P 



^2^ySF{k+^)nk,q)SF{k-^)V{p-k). 

(20) 



Introducing the scalar function 

(i + 7o)r(p,,)ii_M . ^l±p^i^^l^mE), 

(21) 



we write 



{27Tf E + ie-kym^^'^'-^^m^- ^^^^ 



It is easy to see that f (fc, E) satisfies Eq. (fT^ and there- 
fore we can write 



a{E) 



1 



2N, 1 



77 m' 



ImV • V'G"^(r,r',S + ie)|-^p^o- (23) 



Thus to calculate the spectral function in the scalar chan- 
nel, we have to calculate the derivatives of the nonrela- 
tivistic Green's function with respect of f and r' and take 
then the limit f—r' — O, according to ([23]). 



C. Numerical Solution of the Schrodinger Equation 

To obtain the spectral function in the nonrelativistic 
limit, we have to solve Eq. (fT3|) for nonzero e, i.e. com- 
plex energy. We use the numerical method developed in 
Ref. [39] for this purpose which we have extended to in- 
corporate the scalar channel, as discussed in Appendix A. 
At finite temperature all particles have a thermal width, 
but for heavy quarks this is expected to be small. There- 
fore in our study we aim to get the Green's function 
in the limit e 0. In the numerical analysis we used 
Ec = 0.03mc, O.Olmc and O.OOSmc for charmonium, and 
Efc = 0.009TOb, 0.003m6 and O.OOlSm;, for bottomonium. 
For the continuum part of the spectral function all three 
values of the width e give the same result. In the low en- 
ergy part the shape of the spectral function agrees quite 
well for the smallest two ec,6 values. In what follows, we 
will show spectral functions calculated for £c = O.OOBuic 
and ef, = O.OOlSmf,. We note that for bound states close 
to the threshold even a tiny width could have a significant 
effect, namely it could eliminate the bound state peak in 
the spectral function. But on the level of correlators this 
introduces at most a 1% effect. 

For the numerical analysis we need to specify the po- 
tential in Eq. (I13p . The Cornell parameterization of the 



potential turned out to be very successful for the phe- 
nomenological description of the quarkonium spectra, as 
well as a fit Ansatz for the lattice data on quark antiquark 
potential. To include medium effects at high tempera- 
tures, as well as many-body effects at zero temperature 
(e.g. threshold for open charm or beauty production, 
and quarkonia plus glueball production) we will use the 
following parameterization of the potential 



V{r) = 



ar, 



(24) 

The parameters a and a are fixed by zero temperature 
lattice QCD calculations, while other parameters may be 
temperature dependent. We discuss the choice of these 
parameters in the following sections. The potential used 
in the numerical analysis is of course smooth. We used 
a Fermi Dirac function to interpolate between the two 
forms in Eq. (j24p at r = Vmed- At finite temperature 
we also used a more complicated interpolation between 
the short and long distance behavior (see Section ITVl and 
Appendix C). Above deconfinement the singlet free and 
internal energies of static quark antiquark pair have also 
been used as a potential. 



III. QUARKONIUM CORRELATORS AT ZERO 
TEMPERATURE 

In this section we discuss quarkonium correlators ob- 
tained from the spectral functions, which are calculated 
using a potential model matched onto the perturbative 
QCD results at higher energies. 

From the spectral function, determined as discussed in 
section ini we calculate the Euclidean correlators defined 
by 



G(t,T) 



(hj<j{(jj, T)K{uj, T, T). 



(25) 



At zero temperature, K{lj,t,T) = exp(— wr). We com- 
pare the calculated correlators to recent numerical results 
calculated from isotropic [25::] and anisotropic [30] lattice 
formulations. 



A. Numerical analysis of the spectral functions in 
the nonrelativistic limit 

Using Eqs. and we calculate the pseudoscalar 
and scalar spectral functions in the nonrelativistic ap- 
proximation. To do this, we have to specify the parame- 
ters of the potential a, cr, Vmed^ Q?', a' and Vq, as well 
as the charm and bottom quark masses TOc,fc- In this 
work we are interested in QCD with only heavy quarks 
because most of the calculations are performed in the 
quenched approximation. In quenched QCD (pure SU(3) 
gauge theory) the static quark antiquark potential is well 
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m 


Charmonia 
c = 1.19 GeV 




Bottomonia 
mb = 4.575 GeV 


State 


Model 


Lattice 




State 


Model 


Lattice 


I'Po 
2='Po 


3030 
3030 
3437 
3675 
3966 


3012(1) 
3084(1) 
3408(9) 
3739(46) 
4008(122) 




l'5o 
l^Si 
l^Pi 
2^5'i 
2^ Pi 


9406 
9406 
9736 
9874 
10100 


9400 
9426(4) 
9800(16) 
9938(21) 
10181(64) 



TABLE L Charmonium and bottomonium masses in MeV 
calculated in our model and in quenched lattice simulations. 
The values of the charm and bottom quark masses in our 
model are also shown. 



known. In particular, lattice calculations of the potential 
have been extrapolated to the continuum hmit (4^ . It 
turns out that for distances r > 0.4 fm an excellent de- 
scription of the potential calculated on the lattice can be 
given by the Cornell parameterization with a = 7r/12 and 
Here tq is the Sommer parameter 



cr = (1.65-7r/12)r(72 
defined as 



1.65. 



(26) 



As is done in most of the quenched QCD studies, we 
use the phenomenological value of the Sommer param- 
eter To = 0.5 fm. The Cornell parameterization with 
the above parameters gives a fairly good description of 
the lattice data, even at short distances down to 0.1 fm. 
Only at distances r < 0.1 fm the effect of the running 
coupling appear to be important Q . Therefore the Cor- 
nell parameterization is appropriate for describing the 
quarkonium spectrum, which is sensitive to the poten- 
tial in the region O.lfm < r < Ifm. The charm and 
bottom quark masses are chosen such that the poten- 
tial model reproduces the quenched lattice data on char- 
monium [S| and bottomonium spectra j45j . Since the 
NRQCD calculation of Rcf. [i^ does not give the abso- 
lute value of the bottomonium masses, we require that 
the mass of J]b is equal to 9.4 GeV. In table HI we show 
the masses of different quarkonium states and the values 
of the quark masses. In Fig. [1] we show the quarkonium 
spectral function in the pseudoscalar channel for three 
different sets of the parameters /i and rmed- The poten- 
tial ([M)) corresponding to these parameter sets is shown 
in Fig. [2j Here we use a' — a and a' = a. In QCD with 
only heavy quarks, string breaking does not occur until 
distances where the potential becomes comparable with 
twice the heavy quark mass. However, as the energy uj 
is increased the distance between bound state peaks be- 
comes very small (as can already be seen in Fig. [T]). 
Furthermore, at even larger energy, it is possible to cre- 
ate quarkonium plus glueball states which subsequently 
decay into quarkonium states. In this energy regime, it is 
impossible to discriminate between individual states and 
thus the spectral function will have a continuum. This 
will happen at energies of about w ~ 4 — 5 GeV for char- 
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FIG. 1: The nonrelativistic pseudoscalar spectral functions 
calculated for charmonium (top) and bottomonium (bottom) 
using the screened Cornell potential. 



monium and 11 GeV for bottomonium. In this energy 
region the potential model, strictly speaking will break 
down, but the effect of the interaction will still be im- 
portant. To mimic the continuum part of the spectral 
functions, we will choose fj, and rmed such that above 
energies 4 — 5 GeV for charmonium and 11 GeV for bot- 
tomonium, the corresponding spectral functions have a 
continuum. As we will see in the next section, the corre- 
lation function is not very sensitive to the exact choice of 
/i and rmed as long as the continuum threshold is larger 
than 4 GeV for charmonium and 11 GeV for bottomo- 
nium, as shown in Fig. [T] This figure illustrates that in 
the range studied, since the potential is only modified at 
distances larger than 1.5 fm, the lowest lying states are 
not affected by the choice of the parameters rmed and /i. 
In summary, many body effects in the spectral function 
can be simulated by the screened Cornell potential given 
by Eq. (|24p with appropriately chosen /i and rmed- As it 
is discussed in the next subsection this somewhat ad-hoc 
treatment of many body effects has almost no effect on 
the correlation functions in Euclidean time. The width 
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Cornell 

rmgj=2.0fm, |i=0.1GeV 
r^g^=^.5im, |i=0.1GeV 
r^gji=1 .5fm, |i=0.2GeV 



^ 2 I- 

> 
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for tq in the gauge coupling f3 given in Ref. [44| to set the 
scale. As a consequence the value of the lattice spacing 
is smaller; we get a = 0.017 fm for the lattice spacing at 
P = 7.192. Because of this, the meson masses are larger 
than those quoted in Ref. [l^, where the value of the 
string tension ^ = 425 MeV was used to set the scale. 

The relation between the spectral function and the 
nonrelativistic Green's function discussed in Section II 
holds only at leading order. It will be modified by ra- 
diative and relativistic corrections. The radiative correc- 
tions, in particular, turn out to be quite large [47l. lisl. [i^ . 
To take into account these corrections we introduce K 




FIG. 2: The potential at T = for different parameters. 



of the quarkonium states in the spectral functions shown 
in figure [1] are a numerical artifact due to the nonzero 
value of the parameter ec,6. As mentioned in the pre- 
vious section, however, this parameter does not have a 
visible effect on the correlation function. 



B. Direct comparison of the correlation functions 
with lattice data 
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The main focus of this paper is to study the tempera- 
ture dependence of the correlation function. For this the 
form of the zero temperature correlator is not crucial. 
To ensure, however, that the comparison of the lattice 
data and the potential model is meaningful, it is desir- 
able to show that the correlation functions calculated in 
the model agree at least semiquantitatively with the lat- 
tice data. In this subsection we compare the corelators 
from the model calculations to the lattice data at zero 
temperature. Quarkonium correlators have been stud- 
ied in isotropic and anisotropic lattice formulations. The 
correlators of the meson currents calculated on the lat- 
tice require renormalization. The corresponding renor- 
malization constants have been calculated for isotropic 
lattices only: see discussion in Refs. [l^ll^, [M]. There- 
fore, for comparison with our model predictions, we use 
the data obtained on isotropic lattices [1^. We use the 
value of the renormalization constants given in Ref. [1^ . 
In the case of charmonium, we compare our calculations 
to the new lattice calculation on a 48^ x 64 lattice at 
(3 = Q/g^ = 7.192 [46]. In the bottomonium case, we 
compare our calculations to the lattice data of Ref. [1^ . 
In these studies the quark masses, and thus the meson 
masses, were larger than their physical value. For this 
reason, we repeat the analysis from the previous section 
using larger quark masses. In Table |TT] we give the re- 
sulting quarkonium masses as well as the corresponding 
charm and bottom quark mass. In the present analysis, 
we use To = 0.5 fm as well as the interpolation formula 



FIG. 3: The pseudoscalar charmonium correlator calculated 
in our model and compared to lattice data of Ref. [2^. In 
the inset, the corresponding spectral functions a{uj)/uj^ are 
shown. 

factors. These are determined such that the large r be- 
havior of the correlators calculated in our model matches 
the lattice data, i.e. we assume 

2/V 

a{Lj)=K -Im G''' {0,0, E + ie) (27) 

TT 

for the 5'-wave quarkonia and similar relation between 
the derivative of the Green's function for the P-wave 
quarkonia. The values of K are given in Table U for each 
channel. For the study of the temperature dependence 
of the correlator, which is the main objective of this pa- 
per, omitting K would have an effect on the ratio G/Grec 
which is smaller than 0.1%. 

It is important to be aware that the nonrelativistic ap- 
proximation breaks down at large enough energies. On 
the other hand, for large energy, the perturbative result 
for the spectral function given in Eq. ([4]) becomes reli- 
able. Therefore the nonrelativistic spectral function in 
Eq. ([?7|) has been smoothly matched onto the pertur- 
bative relativistic form given by Eq. ^ and as — 0.2. 
In Fig. [3] we show the pseudoscalar charmonium correla- 
tors and the corresponding spectral functions calculated 
for different values of Tmed and ^. The model calcula- 
tion of the correlators show reasonable agreement with 
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2.0 
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11790 
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1.38 
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2.0 




l^Po 


12120 


12295(324) 


2.7 



TABLE II: Charmonium and bottomonium masses in MeV 
from isotropic lattice calculations compared with the masses 
calculated in our model. Also shown are the /^-factors for the 
nonrelativistic spectral functions. 



the lattice data, and the correlator does not depend on 
the values of r^ed and /j,. Furthermore, our calculations 
indicate that even if we use the unscreened Cornell po- 
tential, this causes a change in the correlation function of 
less than 0.5%. Therefore the ad- hoc treatment of many 
body effects discussed in the previous subsection appears 
to have almost no effect on the Euclidean correlators. 

We have used several procedures to match smoothly 
the nonrelativistic function to the relativistic behavior at 
large energies, and have found that the differences in the 
Euclidean time correlator introduced by different proce- 
dures are less than 0.5%. 

To obtain the correct t dependence of the correlators 
at short distances, the relativistic form of the spectral 
function must be used. The nonrelativistic continuum 
leads to a correlator which is more than an order of mag- 
nitude smaller than the lattice data at short distances, 
see Fig. [3l In Appendix B we discuss the additional 
analysis of the behavior of the correlation function at 
T < O.lfm, showing that it is clearly dominated by the 
relativistic continuum contribution of the spectral func- 
tions. This contradicts the statements made recently in 
the literature that the lattice correlator do not 

carry information about the continuum. 

The analysis discussed above for the pseudoscalar char- 
monium correlator we repeated in the vector and scalar 
channels, as well as for the case of bottomonium. Also 
in these channels a reasonable agreement between the 
lattice data and the model calculations has been found. 

Since the Euclidean correlators fall off rapidly with in- 
reasing t it is difficult to judge the agreement between 
the lattice data and the model calculations only looking 
at Fig. [3l It is better to study the ratio of the correlators 
calculated on the lattice to the correlators calculated in 
our model. This comparison reveals some discrepancies 
between the lattice data and the model. These discrepan- 
cies are due to lattice artifacts (at short Euclidean time 
separation), the limited validity of the nonrelativistic ap- 
proach, and in some cases due to the lack of fine-tuning 
of the K-factor, as discussed in Appendix B where we 
show the details of this analysis as well as our results in 
other channels. 



C. Zero temperature spectral function as reference 

In the previous subsection we have shown that the 
nonrelativistic spectral function scaled with proper fac- 
tors, which take into account the relativistic and radia- 
tive corrections, and matched to the perturbative rela- 
tivistic spectral function can give a fairly good descrip- 
tion of the lattice data obtained at quark masses larger 
than the physical values. We have repeated the proce- 
dure described in the previous section for the physical 
value of the quark masses, namely rric — 1.19 GeV and 
TOfc = 4.575 GeV. In this analysis, we used the X-factors 
listed in Table ITll and Vmed — 1.5 fm and fi — 0.2 GeV. 
The calculated spectral function will serve as a reference 
against which the finite temperature results in the next 
sections will be compared. 



IV. QUARKONIUM CORRELATORS AT 
FINITE TEMPERATURE 



As the temperature is increased, quarkonium spectral 
functions will change. This eventually results in a tem- 
perature dependence of the Euclidean correlator G{t, T). 
However, the temperature dependence of G{t, T) is also 
caused by the temperature dependence of the integration 
kernel in Eq. ()25|) . which at finite temperature has the 
form 



K{lu,t,T) 



cosha;(r- l/(2r)) 
sinh (w/(2T)) 



(28) 



To separate out the trivial temperature dependence due 
to the integration kernel following Ref . [1^ , we calculate 
the reconstructed correlator 



Grec{T, T) 



dLua{uj, T = 0)K{uj, t, T) (29) 



and study the ratios G{T,T)/Grec{T,T). If the spec- 
tral function does not change across the deconfinement 
phase transition this ratio will be unity and independent 
of the temperature. In this Section we discuss the the 
pseudoscalar and scalar channels in detail. In the scalar 
channel, there is a zero-mode contribution above the de- 
confinement temperature, i.e. there is a term propor- 
tional to LuS{uj) in the spectral function jHO]. This con- 
tribution is not present in the derivative G'{t,T) of the 
correlator G(r, T) with respect to r. Therefore, in the 
scalar channel we will consider the ratios of the deriva- 
tives G'(r,T)/G;g^(T,r) instead of G(t, r)/G^ec(T, T). 



A. Correlators in the free case 

At sufficiently high temperatures all quarkonium states 
will melt and the interaction between the heavy quark 
and antiquark will be weak. In this limit, quarkonium 
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spectral functions are well approximated by the leading 
perturbative (free field) expression 

-/^ee(c., T) = (^a + 6^ j tanh -^1 - -J 

(30) 

with So = 2mc,b. To obtain an estimate of the tem- 
perature dependence of the correlators, we calculate it 
using a free as the spectral function and consider the ra- 
tio G/Grec- This should provide some upper bound on 
the size of the temperature dependence of the correlators. 
In Fig. m we show the ratios of the correlators and the 
ratios of their derivatives in the scalar channel. In the 
numerical analysis we have multiplied (J free by the factor 
1 -I- Gas /it used in Eq. ([4]) to mimic the leading pertur- 
bative corrections at large oj. This form of the continuum 
ensures that G/Grec approaches one for small r. As one 
can see from the figure, the correlators corresponding to 
free quark propagation are very different from the zero 
temperature correlators, namely they are considerably 
smaller. The differences are considerably larger for bot- 
tomonium than for charmonium and slightly larger in the 
scalar channel than in the pscudoscalar one. 



1.1 
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FIG. 4: The ratio G/Grec for the pseudoscalar and G' /G'r^c 
for the scalar channels for charmonium and bottomonium. 



B. Screening above deconfinement 

Above the deconfinement temperature, we expect color 
screening to take place. The interaction between the 
heavy quark and the antiquark will be modified. We 
assume, however, that the integral equation (fT2|) for the 
nonrelativistic Green's function still holds at finite tem- 
perature but with some temperature dependent potential 
V{r, T). On the lattice color screening is studied by de- 
terminimg correlation function of a color singlet static 
quark antiquark pair separated by some distance r in 
Euclidean time, evaluated at r = 1/T. This is an expec- 



tation value of two Wilson lines 

Tr{L{0)L^{r)) = cxp(-Fi(r, T)/T) (31) 

and defines the so-called singlet free energy Fi (r, T) 
[Slj . Since this object is not gauge invariant, one has 
to calculate it in a fixed gauge. Most of the studies 
of screening of static quarks use the Coulomb gauge 
m m m m, m, instead of using the Coulomb 
gauge, one can insert a spatial transporter between the 
static quark and antiquark, i.e. calculate cyclic Wilson 
loops. Studies of cyclic Wilson loops at finite temper- 
ature have also been performed and have given results 
very similar to the Coulomb gauge calculations [svj . 

Studying Fi (r, T) as a function of the separation r 
three different regions can be distinguished: the short 
distance region, the intermediate distance region and the 
long distance region. At sufficiently short distances, Fi 
is temperature independent and coincides with the zero 
temperature potential. In this region Fi does not de- 
pend on the choice of the correlation function used to 
calculate it (e.g. Coulomb gauge correlator or cyclic 
Wilson loop). At some distance r > r„ied, the singlet 
free energy becomes temperature dependent and deviates 
from the zero temperature potential. The r-dependence 
of the singlet free energy in this region depends on the 
choice of the correlation function (see discussion in Ref. 
[58|). The value of rmed separating the short and in- 
termediate distance regions depends on the temperature 
as rmediT) = 0.43fm/ (T/Tc) tzll- Finally at large dis- 
tances, rT > (1.0 — 1.25) the singlet free energy is expo- 
nentially screened [ist 

Fiir,T)=F^{T)-^^eM-V^iTr). (32) 
3 r 

This feature is independent of the choice of the oper- 
ator. Moreover, Foo{T) is universal and can be ex- 
tracted from the gauge invariant Polyakov loop corre- 
lator {TrL{0)TrL^ {r)) (s^]. This gives the free energy of 
infinitely separated static quark antiquark pairs Foo{T) 
which is smaller than twice the self energy of a heavy 
quark in the medium Voo {T) . This is because the free en- 
ergy contains a negative entropy contribution —TSoo{T) 
|52l |. In fact this entropy contribution dominates at high 
temperatures, making Fao{T) negative for T > 3Tc (52| . 
We would expect that min(Foo(T), 0) is a lower bound 
on Vao{T). Furthermore, for this reason, at temperatures 
close to Tc the singlet free energy provides a lower bound 
on the potential in the sense that V{r,T) > Fi{r,T). 

We assume that the potential V{r, T) shares the gen- 
eral properties of the singlet free energy discussed above, 
i.e. there is a short, intermediate and long distance re- 
gion. In the short distance region, r < r„ied{T) we as- 
sume that V(r, T) is equal to the zero temperature po- 
tential. At large distances, rT > 1.25 we assume that 
V{r, T) has the form 

V{r,T) = VUT)-\—eM-V^iTr), (33) 
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where ai and cti were determined in [43| ■ In the case of a 
screened Cornell potential T4o = cr/^ (c.f. Eq. P^ ). In 
QCD with light dynamical quarks, string breaking occurs 
at distances of about 1 fm. Therefore, assuming /i ~ 200 
MeV we estimate that Voo ~ cr//i — 1.1 GeV, which 
agrees reasonably well with twice the binding energy of 
heavy-light meson 2Ei,ind = 'iMo^s — '^'nicb- Using this 
analogy and realizing that the role of the effective screen- 
ing mass in our case is taken by l/rjned{T), we estimate 
^40(7") — rjnej_{T)a. This way we have specified the large 
distance behavior of the potential. In the intermediate 
distance range rmed{T) < r < 1.25/T we use a Fermi- 
Dirac function to interpolate between the short and the 
long distance behavior, such that the value of the po- 
tential and its first derivative agree at r = rmed and at 
r = 1.25/T. The potential constructed this way is shown 
in Fig. [5] for several different temperatures together with 
the lattice data on the singlet free energy. In Appendix C 
we give details about the construction of this potential. 
Let us note that the above choice of T4o {T) is also moti- 
vated by recent model analysis of the sing let free energy 
with dimension two gluon condensate [59|. 




FIG. 5; The finite temperature potential used in our analysis 
at several values of the temperature together with the lattice 
data on the singlet free energy from Refs. (43,, ,5^, ^l] 



C. Numerical results with F\(r,T) 

From the above discussion it is clear that the singlet 
free energy provides a lower limit for the screened poten- 
tial. Therefore we have analyzed charmonium and bot- 
tomonium spectral functions using the singlet free energy 
as the potential V{r) in Eq. (fT3|) . The numerical results 
for the charmonium spectral functions are shown in Fig. 
[6] (top) together with the corresponding correlation func- 
tions (inset). As one can see, all charmonium states are 
dissolved already at T = 1.2Tc. This is in agreement 
with earlier calculations, which used Fi (r, T) as a poten- 



tial [l6[. The dramatic changes in the spectral functions 
are not reflected in the correlation function, which shows 
only about a 4% enhancement. This is due to the fact 
that even in absence of bound states the spectral func- 
tion is enhanced in the threshold region. The enhance- 
ment near the threshold occurs because the interaction 
between the quark and antiquark is still important at 
this temperatures. An even if the potential is assumed 
to be equal to its lower limit the correlation function does 
not decrease as expected in the noninteracting case. In 
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FIG. 6: The charmonium (top) and bottomonium (bottom) 
spectral functions at different temperatures calculated using 
F\ (r, T) as a potential. The insets show the corresponding 
ratio G/Grcc tog ether with the results from anisotropic lattice 
calculations [33|. G/Grec for charmonium lattice data are 
shown at two temperatures T — 1.2Tc (open squares) and 
l.STc (filled squares), while for bottomonium lattice data at 
l.STc are shown. 

Fig. [6] (bottom) we also show the bottomonium spectral 
function and the corresponding correlators. The ground 
state bottomonium survives as a resonance up to tem- 
peratures as high as 1.5Tc. However, due to the shift in 
the peak position, the correlation functions at this tem- 
perature gets enhanced, resulting in a small increase in 
the ratio G/Grec- The observed enhancement of G/Grec 
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is clearly incompatible with lattice data |28|, |30| (see Fig. 
[S]) . The analysis of Ref . ^6^] which uses the free energy 
as a potential the melting of charmonium IS" states at 
1.2T,. 

We have also calculated the spectral function and the 
correlation functions in the scalar channel. The scalar 
spectral functions show no resonancelike structure above 
the deconfinement temperatures meaning that all P-wave 
quarkonia are dissolved. The interactions between the 
heavy quark and antiquark lead to large enhancement of 
the spectral functions, similar to the one observed in the 
pseudoscalar channel. In Fig. [7] we show the ratio of the 
derivatives with respect to r of G{t,T) and Grecir) for 
charmonium and bottomonium at T = l.5Tc. Similar re- 
sults have been obtained at other temperatures. Our cal- 
culations agree quite well with lattice results within their 
rather large statistical errors. We see again that the large 
enhancement of the spectral function near the thresh- 
old compensates for the dissolution of quarkonium states, 
leaving the correlation function almost unchanged. Thus, 
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FIG. 7: The ratio G' {t) / G'^^^It) in the scalar channel for 
charmonium and bottomonium at T = 1.5Tc. Also shown are 
the lattice data for this ratio from Ref. [30l |. 

contrary to statements made in Refs. [1^, [s^], the dis- 
solution of the IP quarkonium states does not lead to a 
large change in the correlation functions. The weak tem- 
perature dependence of G' [t) / G'^^^{t) was first pointed 
out in Ref. 50]. As we will see in section |Vl the 
larg e change in the P-wave correlators observed in Refs. 
1251 l30l is due to the zero-mode contribution. 



D. Numerical results with the potential V{r,T) 

In this subsection we discuss numerical results from 
our analysis using the potential V{r, T) discussed in sub- 
section IIVBI First, we present the pseudoscalar chan- 
nel. The charmonium and bottomonium spectral func- 
tions at different temperatures are shown in Fig. [5] (top 



and bottom, respectively). Charmonium spectral func- 
tions look similar to the ones discussed in the previous 
subsection. There are no resonance like structures in 
the spectral function; only a large threshold enhance- 
ment. This seems to contradict the conclusions made 
in Refs. [H, H^, [2^, where the analysis of the spectral 
function using the MEM indicated that the first peak sur- 
vives. A more detailed analysis of the spectral functions 
in Ref. [s^l resulted in the more modest conclusion that 
within numerical accuracy no significant temperature de- 
pendence in the pseudoscalar charmonium spectral func- 
tions can be seen. While charmonium spectral functions 
can be reliably reconstructed at zero temperature using 
the MEM, this becomes more difficult at finite tempera- 
ture due to the fact that the extent of the Euclidean time 
is limited to l/T [30]. 

To understand the situation better, in Fig. [8] we show 
our results with the lattice charmonium spectral func- 
tion of Ref. i3Qi] (solid black curve). The first bump in 
the spectral function calculated on the lattice is fairly 
broad and therefore its interpretation as the IS charmo- 
nium state is not obvious. The fact that the bump is 
centered at a; ~ 3.5 GeV instead of the expected w ~ 3 
GeV is a systematic effect. It has been observed that 
also at zero temperature, the position of the first peak is 
shifted toward larger uj when the extent of the Euclidean 
time used in the analysis is limited to small values of 
about Tmax = 0.25 — 0.3 fm [s^]- Lattice calculations 
show, however, that the area under the bump does not 
change within the statistical errors above the deconfine- 
ment temperature. More precisely, the spectral function 
integrated from 2.7GeV to 4.5GeV does not change be- 
tween T = and 1.5Tc. We also calculated in our model 
the integrated spectral function in this interval and have 
found a change of about 5%. Note that the spectral func- 
tion calculated on the lattice has large statistical errors. 
Thus, it is likely that given the statistical accuracy of 
existing lattice data on Euclidean correlators the MEM 
cannot distinguish between a threshold enhancement and 
a true resonance like structure in the spectral function at 
finite temperature. 

In the case of the bottomonium we see only the ground 
state above the deconfinement temperatures, all other 
states are dissolved. In comparison with calculations 
discussed in the previous subsection, we see that the 
medium modifications of the first peak are smaller. This 
is because the potential V{r, T) is deeper than the singlet 
free energy (c.f. Fig. E]). 

In the insets in Fig. [5]we also display the ratio G/Grec 
which shows a much better agreement with the lattice 
data compared to the calculations with Fi discussed 
above. The deviations between the lattice data and the 
results of our calculations in the charmonium case for 
G/Grec are less than 2% for temperatures T < l.?>Tc. As 
the temperature is further increased lattice calculations 
show that G/Grec decreases monotonically [l^HO]- At 
T = 3Tc its value is about 0.9. However, we do not see 
such large decrease in our calculations. We expect that 
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this discrepancy is due to the breakdown of the nonrela- 
tivistic approximation. As the temperature is increased 
the charmonia bound states are melted and therefore the 
typical velocity of heavy quarks becomes of the order of 
the thermal velocity ^jT/rric- At T — 2Tc it becomes 
V ~ 0.7. Thus the nonrelativistic approximation breaks 
down and the simple Ansatz for the nonrelativistic spec- 
tral function given by Ea. (j27p is no longer valid. As 
the temperature is increased beyond this point we expect 
that the spectral function should slowly approach the free 
continuum form (|30p and thus the ratio G/Grec should 
decrease (recall Fig. 2]). We would expect that for bot- 
tomonium the nonrelativistic approximation should be 
valid at higher temperatures and therefore the calculated 
ratio G/Grec should agree better with the corresponding 
lattice data. Results in Fig. [S] show that this is indeed 
the case. 



E. Other choices of the potential 

The choice of the potential discussed above is some- 
what ambiguous. Given the values of rmed{T) and Voo (T) 
there are many ways to interpolate smoothly between 
the short and long distance regimes. We tried differ- 
ent interpolations and found that the spectral functions 
do not depend on the method used. Next, one can ask 
how the results depend on the value of Voo (T) . The free 
energy Fac {T) gives a lower bound on this quantity be- 
cause of the negative entropy term, and the internal en- 
ergy UooiT) provides an upper bound. Therefore we use 
the singlet internal energy calculated in Ref. [6l[ as a 
potential. This quantity has often been used as a po- 
tential when discussing quarkonium properties at finite 
temperature [l^, [l^, |2l[ . The details of the calculations 
are discussed in Appendix D. We find that at 1.2Tc the 
15' charmonium and 25* bottomonium states are present 
as resonances in the spectral functions. At higher tem- 
perature, of about T = lATc the only resonant struc- 
ture which is present is the 15* bottomonium state. We 
also see that the properties of the ground state are sig- 
nificantly modified for this choice of the potential. As a 
consequence, the temperature dependence of the charmo- 
nium correlators does not agree with the lattice results. 
We also find significant deviation in the bottomonium 
correlators at 2Tc. Thus the singlet internal energy is 
not a reasonable choice for the potential. 

The value of rmed{T) in the above analysis was chosen 
according to the analysis of the singlet free energy. There 
are many states which contribute to the singlet free en- 
ergy. As such, the onset of a temperature dependence in 
the free energy at some distance does not necessarily im- 
ply a strong temperature dependence for the potential. 
However, rmed{T) should be smaller than the distance 
where we see exponential screening. Therefore, we take 
rmed{T) = 1.25/T as the upper bound on rmed{T). This 
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The numerical results obtained from the potential con- 



FIG. 8: The charmonium (top) and bottomonium (bottom) 
spectral functions at different temperatures calculated using 
Vi (r, T) as the potential. For charmonium we also show 
the spectral functions from lattice QCD obtained from the 
MEM at l.STc. The error-bars on the lattice spectral func- 
tion correspond to the statistical error of the spectral func- 
tion integrated in the tj-interval corresponding to the hor- 
izontal error-bars. The insets show the corresponding ra- 
tio G/Grec together with the results from anisotropic lat- 
tice calculations [s^l. For charmonium, lattice calculations 
of G/Grec are shown for T = 1.2Tc (squares), 1.5Tc (cir- 
cles), and 2.0Tc (triangles). For bottomonium lattice data 
are shown for T = 1.5Tc (circles) and l.STc (triangles). 



structed with these values of rmed are discussed in Ap- 
pendix D. We see in particular that with such a choice 
for the potential the peaks corresponding to IS" char- 
monium, IP bottomonium, and 2S bottomonium reso- 
nances are present and unchanged in the spectral func- 
tion. The binding energy of these states, i.e. the distance 
between the resonance peak and the continuum, are quite 
small. Therefore thermal fluctuations can destroy these 
states. At l.STc we do not see any resonance like struc- 
tures, except the 15 bottomonium state. Therefore we 
see that for all choices of the potential which are con- 
sistent with the information on color screening coming 
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from lattice QCD, all quarkonium states, except the IS* 
bottomonium, are dissociated at temperatures equal or 
smaller than 1.5Tc. 

In addition, we have studied quarkonium spectral func- 
tions using potentials which are different from the ones 
mentioned above. A common feature of all the potentials 
is that they have the same short distance behavior. From 
the discussion in the previous two subsections it is clear 
that the strong threshold enhancement is due to the short 
distance behavior of the potential. The long distance 
part however, is different. Here we do not use lattice 
data to constrain the long distance behavior of the po- 
tential. The details of the analysis are given in Appendix 
D. The quarkonium correlators obtained with these po- 
tentials also show a weak temperature dependence. More 
precisely, we find that G(t) j Grec(T) in the pseudoscalar 
channel and G' {t) / G'^^^{t) in the scalar channel are tem- 
perature independent and are close to unity. Thus the 
temperature dependence of the correlators does not de- 
pend strongly on the details of the potential. 

V. ZERO-MODE CONTRIBUTION TO 
QUARKONIUM CORRELATORS 

So far, when addressing the scalar channel we have 
discussed only the derivatives of quarkonium correlators. 
The reason for this is the presence of a zero-mode contri- 
bution at finite temperature, i.e. there is an extra finite 
temperature contribution at w ~ in the quarkonia spec- 
tral function 

a, (a;, T) = T) + XI{T)oj5{u:) . (34) 

Here (7^^^^{uj) is the high energy part of the quarkonium 
spectral function discussed in the previous sections, i.e. 
the one at cj > 2mcfi and i = sc, vc, ax for the scalar vec- 
tor and axial- vector channels respectively. For the vector 
channel this has been discussed in Refs. [s^jH^I- In the 
interacting theory the delta function is smeared and has 
a Lorentzian form with the width rj determined by the 
heavy quark diffusion constant D, i.e. r] = T/M/D [6^. 
For values of D which are not too small the contribution 
of the second term in the above equation to the Euclidean 
correlator is given by a constant G'^°"'(r) = Txf (T). The 
susceptibilities xfC^) have been calculated in the free 
theory in Rcf. [35] and read [zl] 

^f*.= (..^)(-|^),,3, 

where Ep = y^p^ + ^^06 ^^'^ — l/(exp(£'p/T) + 1). 
Adding the constant contribution to the P-wave correla- 
tors calculated in the previous section, we can now cal- 
culate the ratio G/Grec in the scalar and axial-vector 



channels. The result of these calculations is shown in 
Fig. [5] both for charmonium and bottomonium (top and 
bottom respectively). Results from isotropic [2^ and 
anisotropic [l^l lattice calculations are also displayed. 
The agreement between our simplified calculations and 
the lattice data is quite good. Our analysis supports the 
observation made in Ref. [s^l that the large increase 
in the scalar and axial-vector correlators is due to the 
low energy contribution to the corresponding spectral 
functions. One should keep in mind that calculations on 
anisotropic lattices were done at quark masses which are 
somewhat heavier than the physical quark masses (c.f. 
Table IhTB)) . In the bottomonium case we see that there 
is a quantitative disagreement between the isotropic and 
anisotropic lattice calculations. This is likely due to the 
fact that the lattices used in Ref. |30| were too coarse 
for precise determination of the bottomonium suscepti- 
bilities xlcaxC^)- ^ constant contribution to the cor- 
relator G'f^{T) = Txi{T) exists at any nonzero tem- 
perature both in the confined and the deconfined phase. 
In the confined phase the quark number is carried by 
heavy charm and beauty baryons (remember that here 
we consider QCD with only heavy quarks) and thus the 
constant contribution is proportional to exp(— 3toc,{)/2^)- 
This contribution is very small. In the deconfined phase 
at sufficiently high temperatures, quark number is car- 
ried by quarks and the constant contribution goes like 
exp(— TOc,b/r), and is described by Eq. which is 

much larger. The fact that we are able to explain the 
behavior of the scalar and axial-vector correlators using 
the ideal gas expressions for the corresponding suscepti- 
bilities (Eq. implies that already at temperatures 
around l.STc, the deconfined charm and bottom quarks 
carry the quark number. This fact directly supports our 
observation that almost all quarkonium states are disso- 
ciated at this temperature. 



VI. DISCUSSION AND CONCLUSIONS 

In this paper we've shown that lattice data on quarko- 
nium correlators and spectral functions may not neces- 
sarily imply survival of different quarkonium states. We 
analyzed quarkonium spectral functions by solving the 
Schrodinger equation for the nonrelativistic Green's func- 
tion of a heavy quark antiquark pair. The nonrelativis- 
tic Green's function is expected to describe the spectral 
function at energies close to the threshold. The results 
of these calculations have been matched to the perturba- 
tive form of the spectral functions at high energies. Al- 
though this matching is not unambiguous the Euclidean 
correlation functions are not sensitive to the details of the 
matching. We have shown that this very simple approach 
can give a reasonable description of the quarkonium cor- 
relators calculated on the lattice. 

Let us note that there are several considerations which 
lead to the identification of the nonrelativistic Green's 
function with the spectral function. One possibility is 
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FIG. 9: The ratio G/Grec in the scalar and axial-vector chan- 
nel at T = l.STc for charmonium (top) and bottomonium 
(bottom). Lattice calculation on isotropic lattices [l^, [l^l are 
shown as filled symbols. Open symbols refer to results from 
anisotropic lattice calculations of Ref. [30l ]. 

to construct an effective theory, the potential NRQCD, 
where the quark antiquark pair is the only dynamical 
field at zero temperature 0, U, ^ . Attempts to gener- 
alize this approach to finite temperatures were recently 
discussed in Refs. [gl, [63|. Bottomonia spectral func- 
tions have been calculated in this approach using Hard 
Thermal Loop perturbation theory resulting in S-wave 
bottomonia spectral functions which are similar to ours 
[63 |. Another possibility to relate Green's functions to 
the quarkonia spectral function is to construct an inte- 
gral equation for the vertex function T{p, q) as discussed 
in section II while systematically taking into account 
medium effects. At sufficiently high temperatures, this 
can be done using perturbation theory. This certainly 
has to be investigated in the future. 

We calculated quarkonium spectral functions at finite 
temperature using a class of screened potentials based 
on lattice calculations of the free energy of static quark- 
antiquark pair. Our analysis shows that independent of 
the details of the choice of the potential, screening effects 
lead to dissociation of all quarkonium states with the ex- 
ception of the IS bottomonium state. In particular, we 
find that 15 charmonium state dissolves at a temper- 
ature below l.STc, contrary to the statements made in 
the literature based on different potential model anal- 
ysis d, [H, iH, [13] as well as on the analysis 
of lattice spectral functions [H [H, HI, HI, [13, IH, . 
We have analyzed in detail the temperature dependence 
of the corresponding Euclidean correlation functions and 



have shown that with some reasonable choice of the po- 
tential, the temperature dependence of the quarkonium 
correlation functions agrees rather well with the pattern 
observed in lattice QCD calculations, i.e. the correla- 
tion functions show very little temperature dependence. 
Thus lattice data on quarkonium correlators and spec- 
tral functions may not necessarily imply survival of dif- 
ferent quarkonium states. This can be explained by the 
fact that even in the absence of resonances, the spectral 
function is significantly enhanced in the threshold region 
compared to the case of free quark propagation. This ob- 
servation may have interesting consequences for quarko- 
nium phenomenology in heavy ion collisions. It means 
that a strong correlation between the heavy quark and 
the antiquark is maintained at all temperatures. The 
heavy quark antiquark pair is created in hard processes 
early in the course of heavy ion collisions and therefore 
the separation between the quark and antiquark is small 
so that some of them will form correlated pairs. The 
above observation implies that the initial correlation be- 
tween the heavy quark and antiquark can be maintained, 
to some extent, through the entire evolution of the fire- 
ball until hadronization. If so, then a fraction of the 
correlated quark antiquark pairs emerging from the hard 
process can form quarkonium states at hadronization. 

In this paper we discussed the correlation functions of 
local meson operators. In lattice calculations, correlation 
functions of extended meson operators have also been 
considered (23l . [65| . These could provide further infor- 
mation on the fate of quarkonium states at high temper- 
atures. However, it is not straightforward to calculate 
these correlation functions in our approach. Also the 
comparison of the lattice data with model calculation is 
less straightforward, e.g. because of the rcnormalization 
issues. 

Lattice calculations of the free energy of a static quark 
and antiquark pair show very strong screening effects 
[4^. On the other hand, the calculations of the S-wave 
quarkonium correlators show very little temperature de- 
pendence [H, [S^l • This seemed to be puzzling. In [3l|, [s^l 
an attempt to describe quarkonium correlators in poten- 
tial model has been made. In those works the spectral 
functions consisting of bound state peaks and perturba- 
tive continuum has been used, and no agreement with 
lattice data has been found. The same conclusion has 
been obtained also in [23 | when perturbative continuum 
has been used. More recently, studies which treat bound- 
and scattering-states on an equal footing were presented 
[2ll . [2^ . In [l^ such treatment improved the agreement 
with lattice in the pseudoscalar channel. In |21j| no agree- 
ment with lattice has been found, which could be due 
to the fact that the internal energy of a static quark- 
antiquark pair has been used as the potential. To obtain 
agreement with lattice for the P-wave correlators inclu- 
sion of the zero-mode contributions are essential [tI]. 
The temperature-dependence of the P-wave correlator 
gives further evidence that most quarkonium states are 
dissolved above deconfinement. In our present work, for 
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the first time, a quantitative understanding of the tem- 
perature dependence of the quarkonium correlators has 
been obtained within a potential model with screening. 
We were able to explain the temperature dependence of 
the ratio G{T)/Grec{T), despite the fact that most states 
are dissolved. Good agreement between the lattice data 
and the potential model prediction was also reported in 
Ref. ,19]. It is important to understand the differences 
between our analysis and that of Ref. [l^. First, dif- 
ferent quarkonium states were reported to exist in the 
quark-gluon plasma. In particular the dissociation tem- 
perature of the IS* charmonium state was reported to be 
1.62Tc [l^. At this temperature, however, the binding 
energy is 0.2 MeV ! A state with such a small binding 
energy cannot be detected in the quark-gluon plasma. 
Second, the authors assumed that only the ground state 
contribute to the quarkonium correlation function. Re- 
cent analysis of the quarkonium spectral function in lat- 
ticeQCD does not support this assumption 30]. Third, 
in [1^ the ratio G{T)/Grec{T) has been calculated as 



G{T,T)IGrec{r,T) 



K{Mg{T),T,T) 
K{Mq{T^Q),T,T) 



(37) 



where Mq{T) is the quarkonium mass and K{uj,t,T) is 
the finite temperature kernel defined in Eq. The 
above equation, however, does not take into account that 
the contribution of the bound state is proportional to the 
square of the wave function at the origin (or the corre- 
sponding derivative) which decreases with temperature 
if screening is present. Thus, in [l9j the agreement be- 
tween the potential model calculations and lattice results 
is accidental and due to the oversimplified approach to 
the problem. 

The analysis presented here can be extended in differ- 
ent ways. Although the IS bottomonium state seems to 
survive in the deconfincd phase, this docs not necessarily 
mean that direct T(IS') production in heavy ion colli- 
sions is not suppressed. In the quark-gluon plasma, the 
IS" state will have a thermal width due to gluon dissocia- 
tion, which, due to the reduced binding energy, could be 
sizable ^66!] . Clearly, the above analysis can be extended 
to full QCD. Such an extension is very timely, as new 
information on screening of static quarks in the quark 
gluon plasma becomes available from large scale lattice 
QCD simulations at almost physical quark masses [s^. 
Finally, it would be interesting to extend the analysis to 
finite spatial momenta. We address these questions in 
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Appendix A 

In this Appendix we extend the method developed by 
Strassler and Peskin in [39] to find, besides the nonrel- 
ativistic S-wave, also the P-wave Green's function of a 
central potential. This method can be used for any po- 
tential which is less singular than 1/r^. 

First, we decompose the Green's function into spheri- 
cal harmonics: 



G'''''{r,r',E + ie) = 

Er=o EL-z ^'('■-;W y,"^(g, <^)[y™(0', </.')]* , (38) 

where gi{r, r' , E + ie) fulfills the Schroedinger equation 



m{E 



V{r)) 



gi{r,r' ,E + ie) = 
= mS{r - r')(39) 



The general solution of this equation can be written as 

giir,r',E + ie) = Aglir>)g'^ir^), (40) 

where g^ , g> are the solutions to the homogeneous 
equation regular at zero and at infinity, respectively, and 
r< = min(r, r'), r> ~ max(r, r'). The constant A is given 
by the Wronskian W (/i, /s; r) = (/i/^ - /[h) \r, 



1 



m 



(41) 



In order to determine these two solutions we then study 
the homogeneous solutions of Eq. ([39]) in the vicinity of 
the regular singular point r — 0. For any potential V{r) 
less singular than 1/r^ there are two solutions, a regular 
(7g(r) and an irregular g{{r), which behave as 



9oir) = r'+' + 

5iM = i + - 



(42) 
(43) 



Since gQ{r) close to the origin has the largest degree, we 
can find a power series solution of the form 



n=0 



(44) 



From this expression we determine the other, linearly in- 
dependent solution gi (r) using standard techniques [68| : 



9{{r)=gi{r) / dr' 



1 



{my 



(45) 



These two solutions are, in general, divergent at r — > c». 
Since by construction they are linearly independent, we 
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may write the solutions and gl. defined in (|40|) as 
linear combinations of 5q(?') and g{{r). Since g[, must be 
regular in the origin, we chose 



9<{r) = dlir) 



where B is defined as 



= - lim 



9oir) 



(46) 
(47) 



(48) 



In order to determine B we solve the homogeneous ver- 
sion of equation p9p numerically. The initial conditions 
for this are determined using equations and (|^5)l . In 
particular, we determine Eq. (j44p up to the fifth power 
and compute 5o((5), <?i(<5), and their first derivatives at 
S = 0.01GeV"\ 

The relationship between and the Green's function 
(for S-wave), or its derivative (for P-wave) is obtained 
from Eq. ((38|) and is 

hm,.,^/^o ImG"'' (r, P,E + iej 
= -m.ltn^^^,^„l^ Sidl^^ (49) 
= -^lim,^oIni(^+i3o) =-i^Imi?o (50) 
lim^,r'^o ImV • V'G"'^ (r, P,E + ie 
^|^lim,^oIm(45"^ 



r dr 



Thus the problem of obtaining the Green's function 
or its derivative in the origin is now reduced to obtain- 
ing the solutions of the homogeneous equations. Care 
should be taken though, since the solutions g[{r) always 
contain terms that introduce divergences as r — > 0. The 
coefficients of these divergences are determined by the 
coefficients of the expansions (|33]) and (|45p . The first 
two coefficients do not depend on the long range part 
of the potential, while the others depend on the details 
of the potential, but can be calculated analytically for 
a given potential not more divergent than l/r^. Since 
the first few coefficients are real, and the 1 /r"-like diver- 
gent terms do not contribute to the imaginary part of the 
Green's function. The divergence of the type g^ir) Inr in 
case of the S-waves is also real and thus it will not enter 
in For P-waves, however, the coefficient of the Inr 

term in ()51|) has also an imaginary part proportional to 
e. Since the width of the quark can be with good ap- 
proximation be taken to be zero, we can take the limit 
of e ^ 0, removing this way the divergence that would 
otherwise enter the imaginary part. In the numerical 
analysis, where we do have a finite width, we solve this 
problem by explicitly subtracting the divergence, which 
we compute from the behavior of gl{r) near the origin. 
The exact form of such terms depends on the exact form 
of the potential. 



Appendix B 

In this Appendix we give further details on the compar- 
ison of the Euclidean correlators calculated in our model 
to those from isotropic lattices [1^, [H, |4^. Since the 
correlation function decays very rapidly with increasing 
T (see Fig. in Fig. [TO] we show instead the ratio of the 
Euclidean correlation functions calculated on the lattice 
and in our model for both charmonium and bottomonium 
channels. In the charmonium case the model is capable 
of describing the lattice data within about 10% accuracy 
for T > 0.2 fm. This is reasonable, as the validity of 
nonrelativistic approximation is marginal this case. At 
smaller separations lattice data deviate from the model 
prediction by about 30%. This is expected due to the 
lattice artifacts discussed in Ref. [s^. For bottomonium 
we have a similar agreement between the lattice data and 
model calculations which extend to smaller r. The rea- 
son for this is that the relativistic continuum part of the 
spectral function is less important for bottomonium, and 
it becomes visible only at smaller Euclidean time sepa- 
rations. Note, that the agreement between the lattice 
data and model calculations in the pseudoscalar channel 
could be improved by fine-tuning the K-factor. However, 
for r < O.lfm the discrepancy between our model calcula- 
tions and the lattice data is larger than for charmonium. 
This is understood, because the 0{ma) discretization er- 
rors are significantly larger for bottomonium. 

We would like to stress again that at sufficiently small 
T the relativistic continuum part of the spectral func- 
tion is important. To demonstrate this point better we 
have subtracted the contribution of the first and second 
peaks in the lattice data using the spectral function re- 
constructed with the MEM and compared the subtracted 
correlator to the one obtained from the free relativistic 
spectral function. This is shown in Fig. 1111 We find a 
reasonable agreement between the free relativistic corre- 
lators and the lattice data. 



Appendix C 

In this appendix we discuss the details of constructing 
the potential V{r,T). It has the following behavior in 
the short, intermediate and long distance regimes: 

Vo{r) = + (rr, r < vq = Vmed 



92 



+ 55'', "Tmed <r <ri 



^i(r-) = Voo- l^xe-^^^-, r > n 

(52) 

with Voo — V{r — > oo) ~ a ■ rj„^^{T). Furthermore we 
have 



0-43/m 1.25 

ro = and ri = Tred = T/Tc. 



T 



(53) 
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FIG. 10: The ratio of the correlators calculated on isotropic 
lattice for charmonium [25l . (top) and bottomonium [2^ 
(bottom) to the model calculations in the pseudoscalar, vector 
and scalar channels. The data in the scalar and pseudoscalar 
channel have been shifted by =1=0.3 for better visibility. 



The four parameters gi, (?2,53 and 174 we chosen such 
that the value of the function g{r) and of its derivative is 
equal to the value of Vb(r) (Vi(r)) and the corresponding 
derivative at r = Tmed {t = fi)- Their numerical values 
are given in Table IIIII The value of gs is zero at T = 
l.bTc, 2.0Ta, while 35 = 0.00689804 at T = 1.2Tc. Now 
we can write the potential for arbitrary r as 

Vir) = fo{r)Vo{r) + {l-fo{r)).h{r)g{r) 

+ (l-/i(r))yi(r), (54) 

where fiir) = ^,r-r- \ with r.i — ro,ri. We choose 

l+exp(^-^j 

S = 0.001. 



Appendix D 

In this appendix we discuss how our results on quarko- 
nium spectral functions and correlators depend on the 
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FIG. 11: The pseudoscalar charmonium correlator calculated 
on the lattice with subtracted resonance contribution com- 
pared to the free continuum correlator. 



T/T, 


To [fm] 


n [fm] 


Ql 


Ql 


91 


92 


93 


94, 


1.2 


0.358333 


1.042 


1.25 


0.5 


-2.601 


2.982 


0.0585 


0.0989 


1.5 


0.2866 


0.833 


0.9 


0.45 


-17.63 


17.94 


-0.234 


0.113 


2.0 


0.215 


0.625 


0.63 


0.35 


-1.058 


1.292 


0.0813 


0.0903 



TABLE III: The numerical values of parameters which deter- 
mine the potential V{r,T). 



choice of the potential. First, we use the internal energy 
U-\ ( r, t) of static quark antiquark pair calculated in Ref. 
p3l| as the potential. We fit the lattice data with the 
following Ansatz 



y(r,T) 



-fir 



-fj,r 



Vo(l 



-fir 



(55) 



The fit parameters are given in Table IIVI In addition we 
use a potential given by Eq. ([M)) for different values of 
the parameters r„ied, A* and Vb and a' — a, a' — a. We 
have analyzed five different sets of parameters labeled as 
IT VI, since the internal energy is labeled as set I. For set 
II we set /i = /i(r) = 2T as predicted by the lattice QCD 
and rmed = '^■2hfm/{T/Tc), i.e. we identified rmed{T) 
with the distance where screening becomes exponential 
(see section ITV Bp . We used the same r^ed in set III but 
/i(T) = l/rmed{T). In parameter sets IV and V we use 
fi — T, while in the parameter set VI we used fj, = 2T . 
The numerical values of Vmed, ^J■ and Vq for the five differ- 
ent sets are given in Table HVl In Fig. [T^lwe show the dif- 
ferent choices of the potentials, including the free energy 
Fi{r,T) and the potential V{r,T) discussed in section 
IV, the internal energy, and the screened Cornell form for 
the two different parameter sets at temperatures of about 
T = 1.2Tc- Below we discuss our numerical results for the 
internal energy and different screened Cornell potentials. 
The quarkonium spectral functions calculated using the 
internal energy as a potential are shown in Fig. 1131 We 
see significant increase in the mass and the amplitude of 
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TABLE IV: The values of the parameter for different sets of 
the potential used in our analysis. 
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FIG. 12: Different potential used in the analysis at 1.2Tc. 
Triangles correspond to lattice data on internal energy from 
Ref. pil| ]. while open squares correspond to lattice data on 
the free energy [5^ ]. 



the 15 state at the lowest temperature for charnionium 
and all temperatures in the case of bottomonium. The 
15' charmonium states dissolves at temperatures around 
1.4Tc, while the corresponding bottomonium states ex- 
ists up to temperatures of about 2.0Tc. The temperature 
dependence of the correlators is shown in Fig. [TH As one 
can see in the Figure, the ratio G/Grec for charmonium is 
significantly smaller than unity for all temperatures. In 
the bottomonium case G/Grec drops below 0.96 in con- 
trast to lattice data. The agreement between the lattice 
data and potential model calculations is worst if we use 
the internal energy as a potential. The quarkonium 
spectral function for the parameter set II are shown in 
Fig. [151 As mentioned above the value of rmed{T) for 
this parameter set is about its maximal possible value and 
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FIG. 13: Charmonium (top) and bottomonium (bottom) 
spectral functions calculated using the internal energy as a 
potential. 



/i = 2r as suggested by lattice QCD. Therefore quarko- 
nium states can survive to higher temperatures and have 
larger binding energies for this choice of the potential. In 
this way we can obtain some kind of upper bound on the 
dissociation temperatures and binding energies for differ- 
ent states. We can see from Fig. [T5]that 15 charmonium 
state as well as 25 and IP bottomonium states survive 
till temperatures as high as 1.2rc for this choice of the 
potential. However, the corresponding binding energies 
are small and therefore interactions with the medium will 
result in a sizeable thermal width and will lead to the dis- 
sociation of these states at this temperatures. At higher 
temperature, namely T — 1.5Tc we no longer see peaks 
in the spectral functions corresponding to these states. 
Thus charmonium 15 state and bottomonium IP and 25 
states are dissolved at temperatures 1.2Tc < T < l.STc 
for this potential. In Fig. [T3] we also show the ratio 
G/Grec for the pseudoscalar channel and G' /G'rec for the 
scalar channel. One can see that this ratio is close to one 
also for this potential. 

We analyzed the spectral functions for other choices of 
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the potential summarized in Table IV and calculated the 
ratio G/Grec which are shown in Fig. [T6]for pseudoscalar 
charmonium. We see that for all potentials we have stud- 
ied this ratio is close to unity. Similar results have been 
obtained for pseudoscalar bottomonium. In the scalar 
channel we have found the ratio G' jG'^^^ is also close to 
unity. 
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The relativistic quark field 'I' can be written in terms of 
nonrelativistic quark and antiquark fields as 'i' = {'ip,x)- 
The simplest operator corresponding to scalar meson 
then is x^^Dip/m [ilj . 
Since the scale in this paper is set h 
value of rmed{T) determined in Ref. [4 
0.48fm/(r/rc) has been rescaled. 

There is a miss-print in Eq. (19) of Ref. [H], (a^'-l-a^')Ji 
should read (a^' + a^')Ji; we thank Edwin Laermann 
and Gert Aarts for clarifications of this issue. 
Since the present paper has been submitted, the role 
of the zero-mode contribution has been also realized by 
the authors of [S^], who extended their analysis of the 
scalar, vector and axial- vector correlators [60] 
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